Existence of axially symmetric static solutions 
of the Einstein- Vlasov system 

Hakan Andreasson 
Mathematical Sciences 
Chalmers University of Technology 
Goteborg University 
S-41296 Goteborg, Sweden 
email: hand@chalmers.se 

Markus Kunze 
Fakultat fiir Mathematik 
Universitat Duisburg-Essen 
D-45117 Essen, Germany 
email: markus.kunze@uni-due.de 

Gerhard Rein 
Fakultat fiir Mathematik, Physik und Informatik 
Universitat Bayreuth 
D-95440 Bayreuth, Germany 
email: gerhard.rein@uni-bayreuth.de 

June 8, 2010 



Abstract 

We prove the existence of static, asymptotically flat non-vacuum 
spacetimes with axial symmetry where the matter is modeled as a colli- 
sionless gas. The axially symmetric solutions of the resulting Einstein- 
Vlasov system are obtained via the implicit function theorem by per- 
turbing off a suitable spherically symmetric steady state of the Vlasov- 
Poisson system. 



1 Introduction 



The aim of the present investigation is to prove the existence of static, 
asymptotically flat, and axially symmetric solutions of the Einstein- Vlasov 
system. This system describes, in the context of general relativity, the evo- 
lution of an ensemble of particles which interact only via gravity. Examples 
from astrophysics of such ensembles include galaxies or globular clusters 
where the stars play the role of the particles and where collisions among 
these particles are usually sufficiently rare to be neglected. The particle 
distribution is given by a density function / on the tangent bundle TM of 
the spacetime manifold M. We assume that all particles have the same rest 
mass which is normalized to unity. Hence the particle distribution function 
is supported on the mass shell 

PM = {OapP^P^ = —c^ and is future pointing} C TM. 

Here gap denotes the Lorentz metric on the spacetime M and if x" are 
coordinates on M, then p°' denote the corresponding canonical momentum 
coordinates; Greek indices always run from to 3, and we have a specific 
reason for making the dependence on the speed of light c explicit. We assume 
that the coordinates are chosen such that 

ds^ = c^goodt"^ + gabdx"'dx^ 

where Latin indices run from 1 to 3 and t = should be thought of as a 
timelike coordinate. On the mass shell can be expressed by the remaining 
coordinates, 

and / = f{t,x"',p^) > 0. The Einstein- Vlasov system now consists of the 
Einstein field equations 

Gap = STTC-^Tap (1.1) 

coupled to the Vlasov equation 

P'^dtf + P'^d^af - V^^yp^dpaf = (1.2) 

via the following definition of the energy momentum tensor: 
rr \ a/2 f f dp^dp'^dp^ 

Tap = c\g\ ' / PaPpJ • (1-3) 

J -Po 

Here \g\ denotes the modulus of the determinant of the metric, and F^^ are 
the Christoffel symbols induced by the metric. We note that the character- 
istic system of the Vlasov equation ()1.2p are the geodesic equations written 
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as a first order system on the mass shell PM which is invariant under the 
geodesic flow. For more background on the Einstein- Vlasov equation we 
refer to [3]. 

In [15l [T71 [18] the existence of a broad variety of static, asymptotically 
flat solutions of this system has been established, all of which share the 
restriction that they are spherically symmetric. The purpose of the present 
investigation is to remove this restriction and prove the existence of static, 
asymptotically flat solutions to the Einstein- Vlasov system which are axially 
symmetric but not spherically symmetric. From the applications point of 
view this symmetry assumption is more "realistic" than spherical symmetry, 
and from the mathematics point of view the complexity of the Einstein field 
equations increases drastically if one gives up spherical symmetry. 

We use usual axial coordinates t G M, p G [0, oo[, z G M, G [0, 2tt\ and 
write the metric in the form 

ds^ = -c^^'^'^'dt' + e^^dp" + e^^^dz" + p^B^e-^'^'^'d^^ (1.4) 

for functions 1^,3,^ depending on p and z. The reason for writing u/c^ 
instead of is so that below i/ converges to the Newtonian potential Un in 
the limit c — >• c«. The metric is to be asymptotically flat in the sense that 
the boundary values 

lim v{p,z)= lim p{p,z) = 0, lim B{p,z) = l (1-5) 

|(p,2)|— >00 \(p,z)\^00 |(p,2)|— >00 

are attained at spatial infinity with certain rates which are specified later. 
In addition we need to require the condition that the metric is locally flat 
at the axis of symmetry, i.e., 

z^(0,z)/c2 +/u(0,z) = lnB(0,z), zGM. (1.6) 

We refer to [4] for more information on axially symmetric spacetimes and 
state our main result. 

Theorem 1.1 There exist static solutions of the Einstein- Vlasov system 
\1.S\) . / fi.3|) with c = 1 such that the metric is of the form and 
satisfies the boundary conditions M.5\) . M.6\) . and the spacetime is axially 
symmetric, but not spherically symmetric. 

It should be pointed out that the above form of the metric excludes solutions 
with non-zero total angular momentum. Since the corresponding general- 
ization induces qualitatively new, additional difficulties it is postponed to a 
later investigation. 
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The strategy of the proof of this result is as follows. Due to the symme- 
tries of the metric the following quantities are constant along geodesies: 



E:= -g{d/dt,p'')=c^e^''^'''p'^ 

= C^e'^/^yi + C-2 (e2/^(pl)2 + e2/^(p2)2 + p252e-2^/c2 (^3)2) ^ (1,7) 



E can be thought of as a local or particle energy and L is the angular 
momentum of a particle with respect to the axis of symmetry. Since up to 
regularity issues a distribution function / satisfies the Vlasov equation if 
and only if it is constant along geodesies, any distribution function / which 
depends only on E and L satisfies the Vlasov equation with a metric of the 
above form. Hence we make the ansatz 



and the Vlasov equation (|1.2p holds. Upon insertion of this ansatz into the 
definition (|1.3|) of the energy momentum tensor the latter becomes a func- 
tional Tq^ = Tai3{i', B, fj,) of the yet unknown metric functions VjB,^, and 
we are left with the problem of solving the field Einstein equations (jl.ip with 
this right hand side. We obtain solutions by perturbing off spherically sym- 
metric steady states of the Vlasov- Poisson system via the implicit function 
theorem; the latter system arises as the Newtonian limit of the Einstein- 
Vlasov system. Our main result specifies conditions on the ansatz function 
(f) above such that a two parameter family of axially symmetric solutions 
of the Einstein- Vlasov system passes through the corresponding spherically 
symmetric, Newtonian steady state. The parameter 7 = l/c^ turns on gen- 
eral relativity and the second parameter A turns on the dependence on L 
and hence axial symmetry; notice that L is not invariant under arbitrary 
rotations about the origin, so if / actually depends on L the solution is not 
spherically symmetric. The scaling symmetry of the Einstein- Vlasov system 
can then be used to obtain the desired solutions for the physically correct 
value of c. 

The detailed formulation of our result is stated in the next section to- 
gether with the basic set up of its proof. The remaining sections of the 
paper are then devoted to establishing the various features of the basic set 
up which are needed to apply the implicit function theorem, and to prove 
various properties of the solutions we obtain. 

We conclude this introduction with some further references to the litera- 
ture. The idea of using the implicit function theorem to obtain equilibrium 





/(x^/) = 0(i?,L) 



(1.9) 
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configurations of self-gravitating matter distributions from already known 
solutions can be traced back to L. Lichtenstein who argued the exis- 
tence of axially symmetric, stationary, self-gravitating fluid balls in this way 
im \T2\ . His arguments were put into a rigorous and modern framework in 
[S]. The analogous approach was used in [TE] to obtain axially symmetric 
steady states of the Vlasov-Poisson system, see also [19]. The approach has 
also been used to construct axially symmetric stationary solutions of the 
Einstein equations coupled to a matter model: In [9] matter was described 
as an ideal fluid whereas in [H [2] is was described as a static or a rotating 
elastic body respectively. Besides the different matter model our investiga- 
tion differs from the latter two in that we employ the rather explicit form of 
the metric stated above and a reduced version of the Einstein field equations 
which closely follows [4]. 

2 Set up of the proof 

In what follows we also use the Cartesian coordinates 



which correspond to the axial coordinates /> G [0, oo[, z G M, c/? G [0,27r]; it 
should be noted that tensor indices always refer to the spacetime coordinates 
t,p,z,ip. By abuse of notation we write i'{p-,z) = u{x) etc. In Section [3] we 
collect the relevant information on the relation between regularity properties 
of axially symmetric functions expressed in the variables x G or p G 
[0,oo[, z G M, respectively. 

We introduce two (small) parameters 7 = 1/c^ G [0, oo[ and A G M. In 
order to obtain the correct Newtonian limit below we adjust the ansatz for 
/ as follows. Let 



(x^,x^,x^) = (/9 cos 93, /9 sin 99, z) G 



so that 



For the particle distribution function we make the ansatz 



/(x,t;) = 0(i?-l/7)V(AL). 



(2.1) 



The important point here is that 



E-lh 



eT'^^^^yrTTH^ - 1 1 
7 '2 
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i.e., the limit is the non-relativistic energy of a particle with phase space 
coordinates {x,v) in case v = Un is the Newtonian gravitational potential. 
For 7 = this limit is to replace the argument of (j) in (|2.ip . We now specify 
the conditions on the functions (j) and tp. 
Conditions on cj) and if). 

{(pi) (j) G C^(R) and there exists Eq > Q such that = for f] > Eq and 
0(ry) > for < Eq. 

{(f)2) The ansatz f{x,v) = cp + U{x)) leads to a compactly supported 

steady state of the Vlasov-Poisson system, i.e., there exists a solution 
U = Un G C^(M^) of the semilinear Poisson equation 

AC/ = AiipN = Att j p + dv, U{0) = 0, 

Un{x) = Un{\x\) is spherically symmetric, and the support of pN G 
C2(R3) is the closed hal\BR,^{fd) where Un{Rn) = Eq and UN{r) < Eq 
for < r < Rn, Uuir) > Eq for r > Rn- 

{p3) 6 + 4-Kr^aN{r) > 0, rG[0,oo[, 

where 

ajv(r) := / cP' (hv\^ + UNir)) dv. 

(V') ^ G C°°(]R) is even with = 1 iff L = 0, and ^ > 0. 
For such a steady state 

lim Un{x) = Un{oo) > Eq. 

)-oo 

The normalization condition Un{0) = instead of Un{oo) = is unconven- 
tional from the physics point of view, but it has technical advantages below. 
Examples for ansatz functions p which satisfy {(pi) and ((/>2) are found in 
[5l[18], the most well-known ones being the polytropes 

cP{E) := {Eq - E)X (2.3) 

for 2 < A; < 7/2; here Eq> Q and (•)+ denotes the positive part. In Section[7] 
we show that for this class of ansatz functions also (03) holds. Numerical 
checks indicate that {(p2>) holds for general isotropic steady states of the 
Vlasov-Poisson system. 

We can now give a more detailed formulation of our result. 
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Theorem 2.1 There exists 6 > and a two parameter family 

('^7,A,-B^,A,/^7,a)(7,A)g[0,<5[x]-5,5[ C C^(]R^)^ 

with the following properties: 

(i) {t^Ofi, Bofi, fiQfi) = ([/at, 1,0) where Un is the potential of the Newto- 
nian steady state specified in {4>2). 

(a) If for 7 > a distribution function is defined by Eqn. \2.1\) and a 
Lorentz metric by ( f with c = l/-^/7 then this defines a solution 
of the Einstein- Vlasov system U.l\) . U.^) . U.3\} which satisfies the 
boundary condition il.6\) and is asymptotically flat. For A 7^ this 
solution is not spherically symmetric. 

(Hi) If for 7 = a distribution function is defined by Eqn. \2.1\) . observing 
\2.S\) . this yields a steady state of the Vlasov-Poisson system with grav- 
itational potential z^o,A which is not spherically symmetric for A 7^ 0. 

(iv) In all cases the matter distribution is compactly supported both in phase 
space and in space. 

Remark. 

(a) The smallness restriction to 7 = is undesired because c is, in a 
given set of units, a definite number. However, if (/, i/, /i) is a static 
solution for some choice of c g]0,oo[ then the rescaling 

f{p,z,p^,p^,p^) = c~^f{cp,cz,cp^,cp^,p^), 
i'{p,z) = c^'^u{cp,cz), 
B{p,z) = B{cp,cz), 
fi{p,z) = n{cp,cz) 

yields a solution of the Einstein- Vlasov system with c = 1. The factor 
in the metric (jl.4p is removed by a rescaling of time. 

(b) The smallness restriction to A means that the solutions obtained are 
close to being spherically symmetric. 

(c) The metric does not satisfy the boundary conditions (|1.5p . but 

lim i^{p,z) = 1^00, lim p.{p, z) = -I'oo/c^ , lim B{p,z) = l. 

|(p,2:)|— i-oo |(p,2;)|— >oo \{p,z)\—>-co 

(2.4) 
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However, if we by abuse of notation redefine v = v — v^o and fj, = 
H + Uqo/c^ then the original condition (jl.Sp is restored and the metric 
(jl.4p takes the form 

= -c^e'^'/'^'cidt^ + 4 (e^^dp" + e^^^dz^ + /,2i?2e-2./-^d(^2^ (2.5) 

with constants ci , C2 > which simply amounts to a choice of different 
units of time and space. By general covariance of the Einstein- Vlasov 
system ([TT]) . p^ . ([L3]) the equations stiU hold. 

(d) In view of [16] part (iii) of the theorem does not give new information 
on steady states of the Vlasov- Poisson system and is stated mainly in 
order to understand the obtained two parameter family of states as a 
whole. However, we note that for the Newtonian set-up in [16] axially 
symmetric steady states were obtained as deformations of a spherically 
symmetric one. The present approach differs considerably from this 
and in principle is more direct. 

(e) In the course of the proof of the theorem additional regularity prop- 
erties and specific rates at which the boundary values at infinity are 
approached will emerge. 

In the rest of this section we transform the problem of finding the desired 
solutions into the problem of finding zeros of a suitably defined operator. 
The Newtonian steady state specified in {(j)2) will be a zero of this operator 
for 7 = A = 0, and the implicit function theorem will yield our result. 
In order that the overall course of the argument becomes clear we will go 
through its various steps, postponing the corresponding detailed proofs to 
later sections. 

The Einstein field equations are overdetermined, and we need to identify 
a suitable subset of (combinations of) these equations which, on the one 
hand, suffice to determine ^,3,1^, and which are such that at the end of 
the day all the field equations hold once this reduced system is solved. We 
introduce the auxiliary metric function 

^ = 7Z/ + /i. 

Let A and V denote the Cartesian Laplace and gradient operator respec- 
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tively. Taking suitable combinations of the field equations one finds that 



V5 ^ 

l\u H • Vv = 47r7 

B ' 



^e(2C-47.) j.^^ + Tn + T22 + ^e2«T; 



33 



+ ^ • V5 = 87r725 (Tn + T22) , 



2pB 



dp{p'dpB)-^d,,B + ^^p{(d, 



1 + P^ I + P^dpi = + 27 pdpvd.v. 



(2.6) 
(2.7) 

(2.8) 
(2.9) 



The last two equations arise from p{Gii — G22) = and pGi2 = re- 
spectively; note that due to (12. ip . Tn = T22 and T12 = 0. Because of the 
asymptotic behavior of B and the structure of the left hand side of (|2.7p we 
write 

B = l + h/p. 

Next, we observe that by taking suitable combinations of ()2.8p and (12. 9p 
we obtain equations which contain only dp^ or dzi respectively, and we 
chose the former. In the above equations the terms T^p are functions of the 
unknown quantities z^, /i, ^ = + p. for which we therefore have obtained 
the following reduced system of equations: 

Aj. = 4^ (^>oo + 7^11 + 7^33) B, e, p; 7, A) - ^V{h/p) ■ Vu, (2.10) 



dpph + d,,h = 8Tr-f^pB<^>ii{i^, B, p; 7, A), 



(2.111 



((1 + dphf + {d,hf) dpi = d,h {d,ph + 272 (p + h)dpud,u) 

+ (1 + dph) (^^{dpph - d,,h) + j\p + h) {{dpi^f - (d.i^f)^ . (2.12) 

We supplement this with the boundary condition l\l.fj\) which in terms of the 
new unknowns and since necessarily h{0, z) = 0, reads 

e(0,z) = ln(l + ap/i(0,z)). (2.13) 

It remains to determine precisely the dependence of the functions ^ai3 on the 
unknown quantities h,^. Since the ansatz (|2.ip is even in the momentum 
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variables p^,p'^,p'^ — the fact that -0 is even is needed here — , all the off- 
diagonal elements of the energy-momentum tensor vanish. The computation 
of its non-trivial components uses the new integration variables 



ri = , s = v , 

7 



the abbreviation 



m{ri,B,v,^) = Be "^'^ \ 
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and yields 

$oo(z^, B, e, p; 7, A) = ^^^''^-^^^^T^o (2.14) 

A foo i'm{r],B,u,'y) 

= — e(2«-47'^) / <^(^)(i + ^^)2 / ^p{Xps)dsdr], 



B 



{e^''-l)/7 



— re ^ / </>(r/) / ip{\ps){m {rj, B,^,"/) — s )dsdri, 



^uiu,B,C,p;j,X) = Tu + T22 (2.15) 

J{ei"-\)l^ Jo 

$33(Z^, 5, e, P; 7, A) = ^^33 (2.16) 

= ^ / '?^(^) / V'(Aps)s dsdi^; 

a J(e7^_i)/^ Jo 

we recall that Tn = T22. The reason for keeping B as argument on the right 
hand sides above is that the matter terms are differentiable in this variable, 
but taking a derivative with respect to h would yield an irritating factor 
1/ p. For elements of the function space chosen below h/ p extends smoothly 
to the axis of symmetry p = 0. 

We now define the function spaces in which we will obtain the solutions 
of the system (j2.10p . (|2.1ip . (j2.12p . As noted above we write, by abuse 
of notation, axially symmetric functions as functions of x G or of p > 
0, z G M; regularity properties of axially symmetric functions are considered 
in Section [3l We fix0<a<l/2 and < /3 < 1, and consider the Banach 
spaces 

Xi := |i/ G C^'"(M^) I i^{x) = v{p,z) = v{p,-z) and \\v\\xi < ooj, 

X2 := {/i G C^'"(R2) I h{p,z) = -h{-p,z) = h{p,-z) and \\h\\x2 < 00}, 

A-s := {e G C2'"(Zk) I C(x) = C{p, z) = e(/9, -z) and < 00}, 
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where 

Zr:= {x eR^ \ p < R] 

is the cyhnder of radius R > 0, the latter being defined in ()2.17p below. The 
norms are defined by 



i\\xo, 



lkllc3.«{M3) + ll(l + k|)^+''Vl/||oo, 
||/i|b4,.(M2) + ||(l + |(p,z)|)2V(Vp)l|oo, 

IICIIc2>'^(z_r), 



and 

A- := A-i X X A-g, \\{y.Ki)\\x := \\v\\x, + \\h\\x, + \\i\W 

Here || • ||oo denotes the L°°-norm, functions in C'^''"(M") have by definition 
continuous derivatives up to order k and all the highest order derivatives 
are Holder continuous with exponent q, 

II II lln- II ^ \D- g{x) - D- g{y)\ 



\o\<k |,j=fc-,s/6K",x^y \x-y 



a 



and denotes the derivative corresponding to a multi- index a G Nq. We 
note that ii h £ X2 then B = 1 + h/ p £ C^{R^), d. Lemma [321 Moreover, 
it will be straightforward to extend ^ to once a solution is obtained in 
the above space. 

Now we recall the properties of the Newtonian steady state specified in 
{<j)2). That condition implies that there exists R > Rn > such that 

UNir) > (Eo + C/7v(oo))/2, r > R. (2.17) 

If 

\W - UnWoo < \Eo - Un{oo)\/4 and < 7 < 70, 
with 7o > sufficiently small, depending on £^0 and Un, then 

g7J/(a;) _ 2 

> Eq for ah \x\ > R. 

7 

This implies that there exists some 6 > such that for all {v, /i, ^; 7, A) gU 
the matter terms resulting from (I2.14p - (l2.16p are compactly supported in 
Bfi{0), where 

U :={{u,h,^;^,X) ex x[0,6[x]-S,6[ \ \\{u,h,0 - {Un,0,0)\\x < S}. 
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In addition we require that (5 > is sufficiently small so that for all elements 
in lA it holds that B = 1+h/p > 1/2, and the factor in front of in (j2.12p is 
larger than 1 /2; since h vanishes on the axis of symmetry, h/p is controlled by 
V/i. Now let an element {u, /i, ^; 7, A) £U he given and substitute it into the 
matter terms defined in (|2.14|) - (j2.16|) . With the right hand sides obtained 
in this way the equations (|2.10p - (j2.12p can then be solved, observing the 
boundary condition ()2.13p and the fact that we require h to vanish on the 
axis of symmetry. We define the corresponding solution operators by 



G'i(z^,/i,C;7, A)(x) 



G3{iy,h,(,;j,X){x) :-- 



1 



1 



\x-y\ \y\ 



Mi{y)dy 



r V{h/p){y)-Vu{y) dy 
47r Jks B{y) \x - y\ 



\n\{p - p,z - z)\ pM2{p, z) dpdz, 



In (1 + dph{0, z)) + / g{s, z)ds, 0< p< R. 
Jo 



Here 



Mi(x) := (^>oo + 7^ii +7^33)(j^(x),-B(x),C(a:),p,;7,A), 
M2{p, z) := j^B{x) $n(i/(x), 5(x), e(x), 7, A), 

M2(p, z) = M2(-p, z) for p < and z eR, and 



g:= {il + dphf + id,hf) 



2^-l 



dzh {dzph + 27^ (/9 + h)dpudziy) 



+ (1 + dph) ( -idpph - dzzh) + ^\p + h) {{dpvf - {dzvf) 



.(2.18) 



Finally we define the mapping to which we are going to apply the implicit 
function theorem as 



F:U^X, {u, h, 7, A) ^ {i^, h, - {Gi,G2,Gs){iy, h, C; 7, A). 

The proof of Theorem 12.11 now proceeds in a number of steps. 
Step 1. 

As a first step we need to check that the mapping T is well defined, in 
particular it preserves the various regularity and decay assumptions. This 
is done in Section [H 
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Step 2. 

The next step is to see that 

J"(C/Ar,0,0;0,0) = 0. 

This is due to the fact that for 7 = A = the choice /i = ^ = triviahy 
satisfies (fCTT) . (f232D . while (l2T0]l reduces to 

Ai/ = 47r$oo(z^,l,0;0,0) 

with 

«>oo(i^,l,0;0,0) =47r^ V2(r? - u) dr^ = j^^cf) Q|f P + z.) d^;; 

notice that /i = imphes that B = 1. By (i?5>2), = Un is a solution of 
this equation, and the fact that [/tv G Afi is part of what was shown in the 
previous step. 
Step 3. 

Next we show that J- is continuous, and continuously Frechet differentiable 
with respect to (z/, /i, ^). The fairly technical but straightforward details are 
covered in Section [5l 
Step 4- 

The crucial step is to see that the Frechet derivative 

L := DI'{Un,0, 0;0,0) -.X ^ X 
is one-to-one and onto. Indeed, 

L{6iy, 6h, (50 = {Su - Li{5v) - L2{6h),6h, 6^ - L^{5h)) 

where 

Li{5u){x) := - I ( -: — - — 7 - 7^ ) aN{y)Su{y)dy, 
L,{Sh)ix) := -3- / V{Sh/p)iy) ■ VUN{y) '^^ 



4vr 7iR3 \x - y\ 

L?X5h){x) := dp6h{0,z) + - [ {dpp6h - dzz5h){s, z) ds, < p < R, 

with otv as defined in {(f)3). To see that L is one-to-one let L{5u, 6h, 6(,) = 0. 
Then the second component of this identity implies that 6h = 0, and hence 
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also (5^ = by the third component. It therefore remains to show that 
= is the only solution of the equation = Li{6i'), i.e., of the equation 

A6iy = iiraNSiy, 6i^{0) = (2.19) 

in the space Xi. Under the assumption on a^r stated in (03) this is correct 
and shown in Section [6l It is at this point that our unconventional normal- 
ization condition in (02) together with the shift in the solution operator Gi 
become important; notice that Li(6i'){0) = 0. 

To see that L is onto let {gi, 92,93) G be given. We need to show 
that there exists {5i',6h,5^) G X such that L{6u,6h,6(,) = {91,92,93)- The 
second component of this equation simply says that 5h = 92- Now 5h G X2 
implies that L'i{5h) E ^Ys, cf. Lemma [3T] (b). Hence we set 5^ = 9^ + L^{5h) 
to satisfy the third component of the onto equation, and it remains to show 
that the equation 

5v - Li{5v) = 91 + L2{5h) (2.20) 

has a solution 5v G Xi. Firstly, L2{6h) £ Xi. The assertion therefore follows 
from the fact that Li : A^i — )• Xi is compact, as is shown in Lemma 16.21 

We are now ready to apply the implicit function theorem, cf. [71 
Thm. 15.1], to the mapping F -.U ^ X; strictly speaking we should suitably 
extend to 7 < 0, but this is not essential. We obtain the following result. 

Theorem 2.2 There exists 81,82 G]0, 5[ and a unique, continuous solution 
map 

S: [0,5i[x]-5i,5i[^ Bs,{UN,m(lX 
such that 5(0,0) = ([/at, 0,0) and 

T{S{-f,X);-f,X) = Ofor all (7, A) G [0,6i[x] - 8i,8i[. 

The definition of J- implies that for any (7, A) the functions (z^, h, ^) = 
5(7, A) are a solution of the equations ()2.10p - (|2.12p . and if / is defined by 
(|2TT]) then the equations ([2S]), dZT]), (IZT2]) hold with the induced energy 
momentum tensor. We can extend ^ to the whole space using the solution 
operator G3 for all x G M^. Also, the boundary condition (II. 6p on the axis 
of symmetry is satisfied: 

C{0,z) = G3iu,h,C)iO,z)=ln{l + dph{0,z))=lnB{0,z); 

recall that ^ = 72^ + /i. For 7 = we conclude first that h = 0, cf. ()2.1ip or 
the G2-part of the solution operator respectively, then the Gs-part implies 
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that ^ = so that the solution reduces to (i^, 0, 0) where v solves 

Az^ = 47r$oo(z^, l,0,p;0,A). 

Since 

$oo(z^, 1,0,/>;0,A) = 47r / </.(r/) / il^iXps) ds di] 

Ju Jo 

coincides with the spatial density induced by the ansatz (12. ip for the New- 
tonian case, cf. |16[ Lemma 2.1], part (iii) of Theorem 12.11 is established. 
If A 7^ then condition implies that / really depends on the angular 
momentum variable L which is not invariant under all rotations about the 
origin, but only invariant under rotations about the axis p = 0. Moreover, 
if the metric were spherically symmetric then the explicit dependence of the 
quantities ^jj on p would imply that the induced energy momentum tensor 
would not be spherically symmetric which is a contradiction. Hence the 
obtained solutions are not spherically symmetric if A 7^ 0. To complete the 
proof of Theorem 12.11 we must show that indeed all the field equations are 
satisfied by the obtained metric ()1.4p . The corresponding argument relies on 
the Bianchi identity VqG"'' = which holds for the Einstein tensor induced 
by any (sufficiently regular) metric, and on the identity VaT'^^ = which is 
a direct consequence of the Vlasov equation (11. 2j) ; Vq denotes the covariant 
derivative corresponding to the metric ()1.4p . The details are carried out in 
Section [8l 

Finally we collect the additional information on the solution which we 
obtain in the course of the proof. 

Proposition 2.3 Let (z^, h, ^) = S{'y, A) be any of the solutions obtained in 
Theorem \2.2\ and define p := ^ — u/c^ and B = 1 + h/ p. Then the limit 
z^oo := liiii|a;|-i-oo i^{x) exists, and for any cr G Nq with \a\ < 1 and x G M'^ the 
following estimates hold: 

|D'^K2;)-z^oo)| <C(l + |x|)-(i+H), 
l)(x)| <C7(l + |x|)~(2+k|)^ 
iD^^ix)] <C7(l + |x|)-(2+IH). 

In particular, the spacetime equipped with the metric is asymptotically 

fiat in the sense that \2.4^ and, after a trivial change of coordinates, also 
ITR) holds. 

Proof. By definition of Gi, limui^^j^ = / dy. The first two 

estimates are established in Lemma [4.2l As to the third one we observe that 
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by the boundary condition (j2.13p and Lemma 14.21 



By (|2.12|) and the known asymptotic behavior of the coefficients in that 
equation which are given in terms of u and h and their derivatives, 

\dpap,z)\< ^ 



l + p^ + \z\ 
cf. Lemma |4.2[ Hence 



1 + \z\^ Jo 1 + s-^ + \z\-^ ~ 1 + \z\^ 

which is the desired estimate for £^{p,z), provided p < \z\. Since we aheady 
know that the metric under consideration satisfies the fuh set of the Einstein 
equations we can now use (I2.8P and (12. 9p to see that also dz^ is given in terms 
of I' and h and their derivatives and satisfies the same decay estimate as dp^. 
Starting from 

\ap,z)\<\c{p,p)\+ l'm{p,s)\ds, 

we can use the decay of dzS, to obtain the decay estimate for ^(p, z) for 
p > z > (or p > —z > 0), and the proof is complete. □ 



3 Regularity of axially symmetric functions 

We call a function / : M'^ — )• M axially symmetric if there exists a function 
/ : [0,oo[xM ^ M such that 

/(x) = /(p, z), where p = \J x\+ x'2 and z = X3 for x G M^. 

In this section we collect some results on the relation between the regularity 
properties of / and those of /. 

Lemma 3.1 Let / : — t- M &e axially symmetric and f{x) = f{p, z) where 
f : [0, oo[xR R. Let k e {1, 2, 3} and a g]0, 1[. 

(a) f £ C^(M^) iff f G C'^([0, oo[xM) and all derivatives of f of order up 
to k which are of odd order in p vanish for p = 0. 
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(h) f is Holder continuous with exponent a G]0, 1[ iff f is. 

Proof. As to part (a) let / E C'^(M^) be axially symmetric. Then / is even 
in xi and X2 and f{p,z) = f{p,0,z). This proves the "only-if" part. For 
the "if" part one checks that the corresponding derivatives of /, which exist 
for /9 / 0, extend continuously to /? = 0. As to part (b) one only needs to 
observe that x i— )• p{x) = ^ x\^ x\ is Lipschitz, since |V/9(x)| = 1. □ 

At several places in our analysis it is convenient to extend functions of 
(p, z) to negative values of p. 

Lemma 3.2 Let h = h{p,z) £ C'^(M?) be odd in p and define 

{ dph{0,z) , p = 0. 

Then b S C^(M'^) and all derivatives of b up to order 3 which are of odd 
order in p vanish for p = 0. By abuse of notation, b G C^(M^). 

Proof. The regularity of b only needs to be checked at p = 0. Since h is 
odd in p it follows that /i(0, z) = dpph{0, z) = for z G M. Hence as p — t- 0, 

b{p, z) = z) - /i(0, z)) ^ dph{{), z), 

and by Taylor expansion, 

dpb{p,z) = -dph{p,z) - \h{p,z) 

= - {dph{0, z) + dpph{T, z)p) 
P 

- ^ (^h{0, z) + dph{0, z)p + ]^dpph{a, z)p^ 

= dpph{T, z) - ^dpph{a, z) (3.1) 
^ ]^dpph{^, z) = 

where cr, r are between and p. All other derivatives can be treated in a 
similar fashion, where one should observe that dzh{0, z) = 0. The regularity 
with respect to x then follows by Lemma |3. 11 □ 
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4 T is well defined 



As a first step we investigate the regularity properties of the functions 
j = 0, . . . ,3, and of the induced matter terms Mi, M2. 

Lemma 4.1 Let (j) and ^ satisfy the conditions ((f)l) and (ip) respectively. 

(a) The functions $00 ^.n-d $33 have derivatives with respect to 1^,5,, P £ 
R and B g]1/2,3/2[ up to order three and these are continuous in 
v,^, B, p,j, X. The same is true for $11 for derivatives up to order 
four. 

(h) For {u,C,h;-/,X) G U, Mi G ^^(M^) and M2 G (72'"(R2) are both 
compactly supported. 

Proof. As to part (a) we note that differentiabihty with respect to and 
p is straight forward. Concerning differentiabihty with respect to and B 
we observe that for j = 0, . . . , 3 the expression ^jj is differentiable once 
with respect to the indicated variables, provided (p G L^^^ proof of 

[16^ Lemma 2.1]. Under the assumption we can first differentiate twice 
before the change to the integration variables r] and s and obtain expressions 
which are essentially of the same form as ^jj, but with or (p" instead of (p 
so that the resulting expression can be differentiated once more. The reason 
why $11 is one order more differentiable is that when differentiating this 
expression with respect to or the integral with respect to s is preserved, 
its integrand is differentiated, and the resulting expression is qualitatively 
of the same type as $00 and can be differentiated three more times. 

Part (b) follows since the functions i^,B,S^ which are now substituted 
into are all at least in C2'"(M^); the fact that ^ is defined only on the 
cylinder does not matter here because the integrals in the definitions of 
^jj yield functions with support in Z/j. □ 

We now show that is well defined, more precisely: 
Lemma 4.2 Let (z^, ^, /i; 7, A) €zU. Then the following holds. 

(a) Gi = Gi(z^, ^, /i; 7, A) G C^'"(M^) is axially symmetric, even in z = X3, 
and\\{l+\x\){Gi-a)\\^, ||(l+|x|)2VGi||oo < 00, wherea = J^dy. 

(b) G2 = G2{v,i,h;^,\) G C^'"(M2) is odd in p, even in z, and 

||(l + |(p,z)|)G2||oc, \\il + \iP,z)\fD^G2\\oo, ||(l + |(p,z)|f I)'G2||oo < 00. 
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Moreover, 



11(1 + \{p,z)\f{G2/p)\\oo, 11(1 + \{p,z)\fD\G2/p)\\o. < oo. 
Here stands for any derivative of order j with respect to (p, z) G M^. 

(c) G3 = G3(z^, ^, /i; 7, A) G C^'"(Zi{) is axially symmetric, eveninz = x^, 
and \\G^\\c2.c.(^Zr) < oo- 

(d) J^{u,(,h;^,X)eX. 

Proof. As to part (a) the potential induced by the matter term Mi, which 
is in C^'°(M^) by Lemma [4. II (b). has the desired regularity and decay prop- 
erties due to standard regularity results in Holder spaces, cf. [131 Thms. 10.2, 
10.3], and the decay of l/|x — y\ and its derivatives together with the com- 
pact support of Ml. As to the source term g = ^{h/p) ■ Vi^ of the second 
term in Gi we notice that i' £ Xi and h £ X2 implies that g £ C^'"(M^) 
with \g{x)\ < C(l + in particular 5 G n L°°(M^). This implies 

the regularity of the potential induces by g and also its decay: 

^'^'^\dy<f ...+ 



^ — y\ •^|^-»/|<|a;|/2 J\x-y\>\x\/2 

<c [ + + A / \9iy)\dy 

J\x~y\<\x\/2 F ~ y\ Fl J 

< C{1 + |x|/2)-3-/' [ + < 1^ 

J\x-y\<\x\/2 \^ ~ y\ fI fI 

for large |x| as desired; for the gradient of the potential induced by g we 
argue completely analogously. 

As to part (b) we first recall that M2 = M2{p,z) is even in p, and the 
actual source term pM2 is odd, compactly supported, and by Lemma [4T] (b) 
and Lemma O (b), M2 G C2'"(M2). Hence G2 G C^'"(M2) is odd in /o G M. 
As to the decay of G2 let suppMa C -Br(O) C M^. Then for \{p,z)\ > 2R 
and (p, z) G supp M2 the estimate 

\ln\{p- p,z - z)\ - ln\{p,z)\\ < 



\iP,z)\ 

holds, and since J pM2 = this implies that 

C 



\G2{p,z) 



G2 [p, z) - 4 / In I (p, z) I pM2 (p, z) dz dp 



- \{p,z)\' 
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the estimates for the derivatives of G2 follow along the same lines. Finally, 
dp{G2/ p) = —G2/p'^ + dpG2/ p which implies that 

C C 

\dp{G2/p){p, z)\ < --^ + -|2|-| . 

\{p^z)\p^ \{p^z)np\ 

This yields the asserted decay when \p\ becomes large. But we can also use 
(j3.1|) to see that \dp{G2/ p){p-, z)\ < G/\z\^. Both estimates together yield 
the asserted decay for dp{G2/p), and the decay for G2/P and dz{G2/p) can 
be dealt with similarly. 

In order to prove part (c) we observe that (|2.18p and the regularity of i' 
and h imply that g and hence G3 E C^'°(Zjj). By construction, dpGs = g. 
Since h is odd in p we find that 

/i(o, z) = dMo, z) = d,Mo, z) = dpphio, z) = 0, 

which implies that 5(0, z) = 0. Thus by Lemma [3Tl G3 G and 
the proof is complete. □ 

5 T IS continuous and continuously different iable 
with respect to z/, /i, ^ 

In this section we give some details of the proof of the following result: 

Lemma 5.1 The mappings 

Gi-.U ^ Xi, i = 1,2,3 

are continuous and continuously Frechet differentiable with respect to v, h, 
and 

Proof. We only show the differentiability assertion and focus on Gi. Defin- 
ing $ = <I>oo + 7*^11 +7*^33 we consider the differentiability only with respect 
to V, and neglecting the dependence on the remaining variables we look at 
the prototype mapping 

G:V^Xi, G{u){x) := [ P^dy, 

Jrs \x - y\ 

where V C is open, <I> € C^(M) and ^ o v has support in a fixed ball for 
all ly gV. Our first claim is that G has the Frechet derivative 

[DG{.)6.]{x) = [ iVMlMy)^^^ . G V, 5. G X,. 
Jr3 \x - y\ 
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In order to prove this claim we need to show that for i/ G V there exists 
e > such that for 5u G Bg(0) C Xi, 

\\G{v + 5v) - G{v) - DG{v)5u\\x, = o{\\Su\\;,^). 

The support property and the standard elliptic estimate imply that 

\\G{u + Su) - G{u) - DG{u)6u\\x, 

<C\\G{u + 6u) G{u) - DG{u)Su\ l^s., 
<C\\^{iy + Si^)- - <^'{u)5iy\\ci. 
< C\\^{i' + 5v) - - $'(z^)(^i^||c2(M3)- 

Clearly, 

W^iy + 5v) - ^{v) - ^'{v)5v\\oo = o(||5i^||oo) < o{\\5v\\x^). 

We need to establish analogous estimates for expressions where we take 
derivatives with respect to x up to second order of the left hand side. Let 
i,j e {1,2,3}. Then 

where both terms on the right are o(||5i^||a'i)- Similarly, 

{Hv + 6v)-^{y)-^'{v)5y) 

= + H - ^"{v) - ^"'{u)5v) d^Md^M 

+ {^'{u + 5v) - <^'{v)) d^^^.Sv + ^"{v + 5u) dx.du d^.Siy, 

and all the terms appearing on the right are odliJi/ll;^!)- This proves the 
differentiability assertion for G. As to the continuity of this derivative. 



\\DG{i^)-DG{D)\\l^;,^^x^)= sup 

IIHUi<i 



/ 



\--y\ 

<C sup ||($'(l/)-$'(f>))(5z/||cl.a(M3) 

IIHUi<i 

< C\\^'{u) - ^'(i>)||c2(R3) as z> z/ in Afi. 



These arguments prove the continuous Frechet differentiability of the first 
part of Gi with respect to v. The derivatives with respect to /i or ^ can be 
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dealt with in exactly the same manner. The source term in the potential 
which represents the second part of Gi can be expanded explicitly in powers 
of 6h and which together with the standard elliptic estimate proves the 
assertion for that term; note that both B and B + 6h/p are bounded away 
from 0. 

The mapping G2 is treated in the same way as our prototype G above, 
except that we have to estimate the source term including its third order 
derivatives, observing that ^>ii has derivatives up to order four with respect 
to B, ^. 

The mapping G3 is easier since the term g defined in (j2.18p can be 
expanded explicitly in powers of Jz^ and 5h where again we observe that the 
denominator in that expression is bounded away from 0. □ 

6 DJ^{Un, 0, 0; 0, 0) is one-to-one and onto 

We recall from Section [2] and Eqn. (j2.19p that in order to prove that the 
map L is one-to-one it remains to show that g = is the only solution of 



in the space Xi. Inspired by the method in [16] we expand g into spherical 
harmonics Yi^, I G Nq, m = —I, I, where we use the notation of [10]; for 
a more mathematical reference on spherical harmonics see [H]. Denote by 
(r, 6, ip) and (s, r, ip) the spherical coordinates of a point x € and y e 
respectively. For I G No and m = —I, . . . ,1 we define 



The symmetry assumptions in the function space Xi imply that gi-i = 
910 = 911 =0) since up to multiplicative constants the spherical harmonics 
with I = 1 are given by sin^e^*"^ and cos^. To proceed, we use the following 
expansion, cf. [TO] . 



Ag = A-KQNg, 5(0) = 0, 



(6.1) 




(6.2) 



1 



00 



47r r[ 



x-y 



EE 



21 + 1 r'+^ 



1=0 m=-l 
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where r< := min (r, s) and r> := max (r, s). In view of ()6.1 



9im{r) = — ^/ / Y——^Y*^{e,ip)dSxaN{s)g{y)dy 
47r /-^ _ rl; 



9/11/ "^(*)-«+T / ^m(^'V')5(y)rf53;'^S 

+ J- Jo r-^ J\y\=s 

= -^TipY aiv(s)-j:j:YS gim{s)ds 

= ~ 21-^1 ( aAf(s)p:pYC/im(s)ds + y aAr(s)^c/im(s) dsj . 
By a straightforward computation we find that gim satisfies the equation 

{r^9L)' = {Kl + 1) +47r r2a^(r)) gi^, (6.3) 

where prime denotes a derivative with respect to r. 

We use this to show that goo = as foUows. We define t(;(r) := 
supo<s<j, |5'oo('S)| so that 1(700(^)1 ^ rw{r); at this point it becomes essential 
that g{Q) = 500(0) = 0. Now (|6.3p can be integrated to yield the Gronwall 
estimate 

w{r) < Att I s\a]^[s)\w{s) ds, r > 0, 



so that w = and hence goo = as desired. 

It therefore remains to consider gi^ with / > 2. For these we prove the 
following auxiliary result. 

Lemma 6.1 Let a € Cc([0, cx)[) and A > 6e such that A + 47rr^a(r) > 
for r G [0,oo[. Let u S C^([0, oo[) he a hounded solution to 

(r^n')' = (A + A'Kr'^a{r))u. (6.4) 

Then u = 0. 

Proof. We fix > such that a(r) = for r > r^. Multiplying ()6.4p with 
u and integrating by parts we obtain for r > 0, 

{\ + A^s^a{s))u^{s)ds= / {s^ u {s))' u{s) ds 

Jo 

/•r 

= r^u'{r)u{r)- I s^{u'{s)fds. (6.5) 



23 



Now if there exists tq > so that u(ro) = or n'(ro) = then (j6.5p imphes 
that n(r) = u'{r) = for r G [0, rg]. The unique solvabihty of (j6.4p for 
r > To then shows that n = as claimed. 

So we assume now that u{r) ^ and n'(r) / for r > 0. Since (|6.4p 
is invariant under u — )• — u, we may suppose that u{r) > and n'(r) > 
for all r g]0,oo[; note that ()6.5p enforces uu' > on ]0, oo[. For r > > 
()6.4p simplifies to {r'^u')' = Xu, which has the solution 

/ N _ (/ + l)u{ra) +rau'{ra) ( rV lujrq) - rau'{ra) ( Va 
(2/ + 1) \rj + (2/ + 1) [r 

Therefore u is unbounded which is a contradiction. □ 

Since g £ Xi, Eqn. (|6.2p implies that gim is bounded. Due to (cpS) we 
can apply Lemma ETT] to conclude that gim = for alH > 2, and thus g = 
as desired. 

We now prove the compactness result which was needed to show that L 
is onto. 

Lemma 6.2 The mapping K : Xi Xi, 

{Kw){x) = [ ^^^M^dy 
Jr3 \x - y\ 

is compact. 

We remark that the operator Li has the form Li{5i'){x) = —K{5v){x) + 
K{5v){Q) and is compact if K is, since the mapping v i— )• z^(0) is continuous 
on Xi. 

Proof. First we observe that the mapping 

[ (p' ( -\vf + u] dv = 2V27r [ (t)'{E) \/E -udE 

is in C^(M), and since Un e C^(M^) the function un is in C'^{R^). Hence 
afyfW S C^'^^'^{M^) for any w S Xi, and since a < 1/2 the mapping K is well 
defined. 

We fix a function x £ C^(M'^) such that < x < Ij x{^) = 1 for |x| < 1, 
and x{x) = fo^^ l^^l ^ 2. Let xni^) = x{^/R) for i? > and define 

{Krw){x) = xr{x){Kw){x). 

We show that Kn — ?■ -fC in the operator norm as i? — t- oo. To this end, let 
Cij = 1 — Xr so that for w £ Xi and x G , 

{Kw - Krw){x) = Cr{x){Kw){x), (6.6) 
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and the latter vanishes for |x| < R. Now let < 1- For cr G Nq with 

\a\ < 3 it follows that 

iD^'iKw - Krw){x)\ < CrXx) \D''{Kw){x)\ 

+ \crD^CRix)D--HKw){x)\ 

0<T<a 

' ' 0<T<a 

constants denoted by C do not depend on x or R. In order to estimate 
the Holder norm of D'^{Kw — Krw) for \a\ = 3 we take x,x £ with 
I a; I > I a; I and again apply the product rule to the expression (|6.6p . Adding 
and subtracting terms we have to estimate expressions like 

\{D\n{x) - D-Cr{x))D''-^Kw){x)\ < ^\x - x\ 

and terms like the following: 

iD^CRix)] [ \D''-^\y\-^ \{aNw){x - y) - {aNw){x - y)\ dy; 



if t| = 3 we throw one derivative onto aNW. The latter quantity together 
with its first order derivatives is Holder continuous. The factor in front of 
the integral vanishes for |x| < R, so we need only consider \x\ > \x\ > R. 
Since the domain of integration extends only over y with |y — 2;| < R]\f or 
|y — x| < Rn we can on the domain of integration estimate \y\ > \x\ — 
\x — y\ > R — Rn ^ -R/2 or analogously with x instead of x, where we 
assume that R > 2Rn. Since |Z)'^~'^|y|^^| < \y\^^ with j > 1 the term 
under consideration can be estimated by CR~^\x — x\'^ and altogether we 
conclude that 

\\Kw - Kjiw\\c3,c,(K3) < C/R. 

Recalling the definition of the norm || • we see that the following chain of 
estimates finally shows that Kn K m. the corresponding operator norm 
as desired: 

\V{Kw){x) - V{Krw){x)\ < Cr{x) [ '"^(^)^yi dy 

7r3 \x-y\^ 



+ R-^Vx{x/R)\f \^-^y)^^)\ dy 
Jr3 \x - y\ 

C _i C 

< C{l + \x\)-^-^R-'~^-^\ 
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To complete the proof we have to show that Kji is compact for any R > 
on the space Xi. First the fact that oat G C|(]R^) impHes that 



is continuous, and the same is true for 

Kr : ^ C^'^/\BsRm 

where we note that all the functions Krw with w S C3'"(M'^) are supported 
in B3r{0). Since a < 1/2 the embedding 

C3'1/2(53h(0)) C3'"(S3k(0)) 

is compact, and because of the support property we conclude that 

Kr:Xi^ Xi 

is compact; on VKrw the weight + only amounts to multiplication 

with a bounded function. □ 

7 Discussion of Condition (03) 

In this section we investigate Condition (03) for the case of the polytropic 
steady states (j2.3p . We first allow for the general range k £] — l/2,7/2[ of 
polytropic exponent. Using the elementary integration formula 

s-hv\^Ydv = i2nf/^^p:^s';K seR, (7.1) 



2' W+ ' ' r(A;+|) 
the Poisson equation in {(f)2) is found to be 

for Un = U]\fir). According to [18] there exists a solution Un such that 
Un{0) < Eo, U'j^iO) = 0, Un{Rn) = Eo, UN{r) > Eq for r > Rn, and 
U'j^{r) > for r e]0, Rn[- For z := Eq — Un this means that 

_J^(rV)' = 47rc„z!^, where n := A; + ^ g]1, 5[, Cn := (2^)^/^ + , 
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and furthermore z(0) > 0, ^'(0) = 0, z{Rn) = 0, and z'{r) < for r G 
]0,i?Ar[. In terms of z the function ajy from (03) reads 

a^ir) = -i2.f^^^^zirf_;'^ = -nc^z{r)r\ 

where once more (|7.ip was used. Thus condition ((/)3) is equivalent to 

^'Kncnr'^z{r)1'^ <Q. (7.2) 

Now consider the function C(s) := z{as) for a := (47rc„)~^/^. It is found to 
satisfy the Emden-Fowler equation 

- \{s\')' = a (7.3) 

and C(0) > 0, C'(0) = 0, C(so) = for sq := Rn/oi, as weh as C'{s) < for 
s g]0, so[. In terms of s = a~^r condition (|7.2|) becomes 

s\{s)l-' < -. (7.4) 
n 

The left-hand side can be conveniently expressed by means of the dynamical 
systems representation of (17. 3p . For, let 

^(t):=-^>0, V{t):=-'-^>0, t:=lns, 

where we consider t g] — oo,lnso[. Then 

if = U{3-U -nV), V = V{U + V -I), (7.5) 

and U{t)V{t) = s^C(s)+~^, which provides the relation to (j7.4p . Thus we 
have to verify that U{t)V{t) < 6/n. In the terminology of [6, p. 501], 
where m = 0, C is an £^-solution to ()7.3p . Thus O Prop. 5.5] implies that 
{U{t),V{t)) lies in the unstable manifold of the fixed point P3 = (3,0) of 
()7.5p . In particular, we have limt^_oo(f^(i), ^(i)) = (3,0). Also note that P3 
is of saddle type with eigenvalues —3 and 2; the corresponding eigenvectors 
are (1, 0) and (— 3n/5, 1). Since the line V = — separates the regions 
U > (below the line) and U < (above the line), a phase plane analysis 
reveals that we must always have U{t) < 3, so that W{t) := U(t)V{t) < 
3V{t). In addition, it is calculated that V and W are solutions to the system 

V = V{V -1) + W, W = W{2-{n-l)V), (7.6) 
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such that hmt_i.„oo(^(i)) W{t)) = (0, 0). The origin is a fixed point of saddle 
type for (j7.6p , the eigenvalues are — 1 and 2 with corresponding eigenvectors 
(1, 0) and (1, 3). Note that 14^ > for F < W < for V > V >0 
above the curve V i— )■ ^(1 — V), and ^ < below this curve. Since the 
curve has unity slope at V = 0, it follows that (y(t),W{t)), lying in the 
unstable manifold of the origin, will be above the curve for t very negative. 
Then a phase plane analysis shows that this property persists for all times. 
In particular, we always have ^ > 0, and W is increasing until it reaches its 
maximal value for to such that V(to) = Thus our original problem of 

proving (03) is equivalent to showing that W^Iq) = maxW < 6/n. Thanks 
to the preceding observations the parametrized curve t i— )• {V{t),W{t)) for 
t G] — oo,to] can be rewritten as a curve W = W{V) in the (F, Ty)-plane 
which solves 

dW W{2 - (n - 1)V) 



dV V{V-l) + W ' 
and which is such that VF(0) = and W{:^^) = maxM^. 



(7.7) 



Lemma 7.1 If k < 7/2 is sufficiently close to 7/2, then ((j)2>) holds for 
given by {H 



Proof. If W{V) < 1 < 6/n for all V g]0, then we are done. Hence 
we assume that W{Vo) = 1 for some Vq G]0, Then 1 = W{Vo) < SVq 

yields Vq > 1/3. Since W{V) > I fox V > Vq, it follows that V{V-l) + W = 
{V -1)^ + V + W -1>V, so that by (17^ . 



In maxT^ = / < / ^—^dV 

Jw=l W - Jy^ V 

< r'^i^^^^^dv 

A/3 V 

= 21nf— -2 + -(n-1). 
Vn-1/ 3^ ^ 



Therefore 



maxW^<- ^exp(-(n-l)-2). (7.^ 

(n — 1)^ \6 



At n = 5 the relation 



9 ,-2/3 . 6 



4 5 

holds. Hence it follows from (j7.8p that maxVF < 6/n is verified for n 
sufficiently close to n = 5. □ 
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The method of proof for the preceding lemma can be refined as follows. 
Fix A<6/n. Then W{V) < A iov V e [0, would be acceptable. Hence 
we can assume that W{Vo) = A for some Vq g]0, Then A = W{Vo) < 

3Vo shows that Vq > A/3. From W{V) > ^ for ^ > Vb we obtain 

ln(max VK) — In A 

.maxH/^^ /-^ (2- (n- ^-^ 

J A/3 V^-V + A 
5 — n 



n \ / 3-2A 

arctan I , + arctan 



V4yl- 1)/ \3y/AA-lJ _ 



/ 6- ZA \ 
\3JAA-l) 



n-l\, f9\An^ -2iA + l)n + Q + A] 
In 



2 J \ A{A + Q){n-l 
Therefore max W < (n) where 



X exp 



9[An2 - 2{A + l)n + 6 + A] 

n ( 5-n \ ( 3-2A 

arctan , + arctan 



1)^ V3V4^ - 1 



For different yl it can be checked (e.g. using Maple) for which values n G 
]1, min{6/^, 5}[ the relation ^ j^{n) < 6/n is verified. Taking j4 = 1 we 
get at least n £ [2.6, 5[, for A = 6/5 we get at least n £ [2.35,4.85], and 
for A = 2 we get at least n G [2.1,2.5]. In summary, the desired relation 
maxM^ < 6/n can be obtained for at least n G [2.1, 5 [, which corresponds to 
at least k G [0.6, 3.5[ in (|2.3p . Notice however that the regularity assumption 
on (j) requires k > 2. 



8 The field equations iiold 

For a metric of the form (II. 4p the components 00, 11, 12, 22, and 33 of the 
field equations are nontrivial. We have so far obtained a solution i^,B,^ of 
the reduced system (12. 6p . (12. 7p . (I2.12P where the appearing components of 
the energy momentum tensor are induced by a phase space density / which 
satisfies the Vlasov equation (II. 2j) . We define Eaj3 ■= Gap — 87rc~^TQ^ so 
that the Einstein field equations become E^jB = 0. By (j2.7p . 

Eii + E22 = ^. (8.1) 
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Using this information (j2.6p says that 



or 



^2g4./c-^00^ ^2^2^33 ^0. (8.2) 

The Vlasov equation imphes that S/aT°^^ = 0, and VoG"^ = due to the 
contracted Bianchi identity where Vq denotes the covariant derivative cor- 
responding to the metric (jl.4p . We want to use these relations to show that 
the remaining components of Ea^ vanish also, but there is a technical catch: 
The metric, more specifically ^, is only C^. To overcome this complication 
we approximate ^ by functions The induced Einstein tensor Gn^ 
again satisfies the Bianchi identity. Taking /3 = 1 and letting n — )■ oo we 
obtain the equation 

d,E^^ + (^4a,/i + ^) E^^ - pBe-^i{B + pd,B) = q, (8.3) 

where (18. 2p has been used to eliminate E^^ and we recall that ^ = u/c^ + jj,. 
Here dzE^'^ is at first a distributional derivative, but since all other terms 
in the equation are continuous this derivative indeed exists in the classical 
sense. The same approximation maneuver can be performed for /3 = 2 to 
obtain the equation 

d,E'' + (4d,p + i + ^) E'' - p'Be-'^d^B E^' = (8.4) 

which holds for p > 0. However, if we multiply this equation with p we 
obtain an equation which holds for p > 0. This is because E^'^{0,z) = 
which is nothing but the boundary condition (j2.13p on the axis of symmetry 
which we have incorporated into our integration of (j2.12p . We eliminate E^^ 
from ()8.3p . (j8.4p and write the resulting equation for E^"^ in terms of 



X := pe^^'BE^^. 

The result is the equation 

B + pdpB 

which again holds for p > 0. Since -'^(O, z) = and since any characteristic 
curve of this equation intersects the axis of symmetry p = we conclude 
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that X vanishes identicahy. By (j8.3p the same is true for E^^ so that -E^^ = 
^33 ^ ^00 ^ Finally we observe that by (|2J2]1 . 

(0 B\ 3 B 

1 + P^j {Ell - E22) + P^Ei2 = 0. 

Since E12 = this means that En = E22, and with (|8.ip we conclude that 
£^11 = £'22 = 0, and all the non-trivial field equations are satisfied. 
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